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Abstract. We show how to create maximally entangled dressed states of a 
weakly interacting multi-partite quantum system by suitably tuning an external, 
periodic driving field. Floquet theory allows us to relate, in a transparent manner, 
the occurrence of entanglement resonances to avoided crossings in the spectrum of 
quasi-energies, tantamount of well-defined conditions for the controlled, resonant 
interaction of particles. We demonstrate the universality of the phenomenon 
for periodically driven, weakly interacting two-level systems, by considering 
different interaction mechanisms and driving profiles. In particular, we show that 
entanglement resonances are a generic feature of driven, multi-partite systems, 
widely independent of the details of the interaction mechanism. Our results are 
therefore particularly relevant for experiments on interacting two-level systems, 
in which the microscopic realization of the inter-particle coupling is unknown. 



PACS numbers: 03.67.Bg, 42.50.Dv, 32.80.Qk, 42. 50. Hz 

1. Introduction 

Strongly entangled multi-partite states of a collection of quantum mechanical objects 
are a central resource in quantum information science, think, e.g., of one-way quantum 
computing [T. Beyond that, entangled many-particle states naturally occur in 
quantum phase transitions [2J, and in the condensed phase [3], and recently they 
are also debated in the context of non-equilibrium processes in biological tissue 0J- 

A viable way to achieve long-lasting entanglement in autonomous quantum 
systems is to adjust the system parameters such that one of the eigenstates becomes 
entangled, and then to prepare the system in this eigenstate [5]. The appeal of this 
scheme is that, in the idealized situation of an isolated quantum system, the system 
remains in this eigenstate in the course of time, and hence preserves its entanglement. 
Even if the idealization of an isolated system is not justified, the coupling to a noisy 
environment, if not too strong, typically affects system eigenstates least [6], rendering 
their entanglement inherently robust against decoherence. Accordingly, entanglement 
of eigenstates of autonomous quantum systems has been studied extensively in the 
literature, e.g., for spin chains [7]. 

In this paper, we extend this scheme to periodically driven quantum systems, 
since the application of external, time-periodic driving fields is often the most 
elementary way to control a quantum system - think of, e.g., trapped ions [5], cold 
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atoms 9 , or color centers in diamonds [10]. In the presence of a driving field, 
the system Hamiltonian is time- dependent, and the notion of eigenstates becomes 
obsolete. However, for time-periodic driving, the dressed state, or Floquet, picture 
applies, constituting a framework that is essentially equivalent to the concept of 
eigenstates of autonomous systems [TTJ [T2 H3J. In particular, it provides quasi- 
stationary solutions of the dynamics, i.e., solutions for which the system dynamics 
repeats itself periodically after each driving cycle. Thus, if one of the Floquet states 
is strongly entangled over the entire driving period, and if the system is initially 
prepared in this state, its entanglement is preserved - in complete analogy to entangled 
eigenstates of autonomous systems. Nevertheless, beyond this formal analogy, Floquet 
states can differ strongly in their properties from the eigenstates of the respective 
undriven system. This is manifest in micromaser physics |14[ I15j . rapid adiabatic 
passage [IB], the coherent destruction of tunneling [T7], non-dispersive wave packets 
[T5] , scenarios of Anderson / dynamical localization in light-matter interaction [18 , 
and the asymptotic persistence of quantum coherence under environment coupling 

In this paper, we show how maximally entangled Floquet states of interacting 
two-level systems can be generated by controlling the parameters of an external, time- 
periodic driving field. Taking a perturbative approach in the interaction strength 
between the qubits, we connect their occurrence to avoided crossings of many- 
particle states in the quasi-energy spectrum. This is conceptually closely related to 
controlled iV-particlc interactions that have recently been demonstrated in a cold 
Rydberg gas [21], only that here we assume alternating instead of static control fields, 
what enormously widens the versatility of the control strategy. For autonomous 
quantum systems, the connection between avoided crossings and entanglement has 
been discussed earlier, as, e.g., in spins chains [55] [23]. Our generalization to 
periodically driven system opens up the possibility of creating strong, stationary 
entanglement in multi-partite quantum systems, by simply tuning the parameters 
of an external driving field. 

2. Theoretical framework 

2.1. Floquet theory 

We start with a recollection of Floquet theory, on which we rely throughout this work, 
and which is the semiclassical [f] variant of dressed state theory [TTJ [T2] . It is founded 
on the theorem that, given a Hamiltonian H(t) = H(t + T) with period T, every 
solution of the Schrodinger equation 

i a !*(*)> = #(*)!*(*)) (s = i) (i) 

can be written as a superposition 

|*(<))= ^e^^i)) (2) 

i 

of mutually orthogonal, T-periodic Floquet states |3>i(£)) = \$>i(t + T)). The quasi- 
energies appearing in the phase factor are real numbers, and the time-independent 

X Floquet theory is semiclassical in the sense that the driving field is not treated in a quantized 
fashion, but as time-dependent term in the Hamiltonian. 
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weighting factors of the superposition are Oj = ($^(0)1^(0)). Hence, once the Floquet 
states and quasi-energies of H{t) are known, the time evolution of the system is 
available for arbitrary times and arbitrary initial states. 

The Floquet picture can be regarded as the generalization of the concept of 
eigenstates to periodically driven systems: Consider, for the moment, an autonomous 
Hamiltonian Ho, plus a weak periodic perturbation V(t) — V(t + T), such that 
H{t) = Hq + XV (t). In the limit A — > 0, a possible choice of Floquet states and 
quasi-energies are, respectively, the eigenvectors |$$) and eigenvalues J5{ of H , which 
can be found by solving the eigenvalue problem Ho = Ei\$i). In this case, 
Eq. ([2| is simply the well-known time evolution of eigenstates, reflecting the fact 
that an autonomous system will not leave an initially prepared eigenstate. For finite 
A, the Floquet states of H(t) become time-dependent. As long as the driving V(t) 
varies only slowly in time, the system dynamics will be adiabatic; i.e., the Floquet 
states are the instantaneous eigenstates of H(t), given by H(t) \<&i(t)) — Ei(t) 
However, for fast (and/or strong) driving, the adiabaticity assumption is not justified 
(cf. Refs. [231 [25] for a detailed discussion) , and Floquet states no longer coincide 
with the instantaneous eigenstates of H(t). Nevertheless, Eq. ([2| still implies that 
after initially preparing a Floquet state \$k(t)) (i.e., a, = So-), the system remains 
in this state for arbitrarily long times, and only gains a dynamical phase e~ l£kt . The 
expectation values of any (time-independent) observable then exhibits at most T- 
periodic time-dependence. 

As in the static case, Floquet states can be transferred into each other via 
frequency-matched probe fields [26], or by adiabatic passage [13]. An important 
difference to the concept of eigenstates of autonomous systems is, however, that 
Floquet states and quasi-energies are never unique: Introducing the driving frequency 
uj = , one finds that with every Floquet state | $<(*)), e~ ikuJt | $;(£)) is an equivalent 
Floquet state, for arbitrary integer k. This is because the latter expression retains 
T-periodicity and, after inserting it into Eq. ^ and shifting the quasi-energy Si by kw, 
leads to the same solution |\&(£)) of the Schrodinger equation [27]. (More formally, 
this equivalence relation defines a rest class structure in the set of Floquet states 
[25].) As a consequence, the spectrum of quasi-energies is identical in every interval 
[kuj, kuj + uS). By restricting quasi-energies to a single such Floquet zone, one gets rid 
of the ambiguity and is left with complete and mutually orthogonal set of d Floquet 
states, d being the dimension of the Hilbert space that H(t) is acting upon. 

2.2. Floquet state entanglement 

In this paper, we study the entanglement of Floquet states of a periodically driven, 
closed quantum system, with the following motivation: If the system is initially 
prepared in one of its Floquet states |$j(0)), the solution of the Schrodinger equation 
is, according to Eq. |*(i)) = e" fei * |$i(t)). Since a global phase factor is irrelevant 
to entanglement [|] the system remains entangled as long as |$t(f)) is. At best, \&i(t)) 
is maximally entangled over the entire driving period t € [0, T). The system then 
remains maximally entangled for, in principle, arbitrarily long times. For this reason, 
it is desirable to understand under which conditions Floquet states are strongly, maybe 

§ This argument also guarantees that all Floquet states of the same rest class have identical 
entanglement properties, since they differ only by a phase factor e ~' kLjt . Hence, entanglement 
properties of Floquet states do not depend on the Floquet zone of the spectrum one is looking 
at. 
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even maximally, entangled. 

Of course, under realistic conditions, decoherence will affect the system after a 
certain time and lead to a deviation from the perfectly coherent scenario described 
by a pure quantum state, that we study in this work. Nevertheless, even in presence 
of weak environment coupling, Floquet states are typically the most robust states of 
a periodically driven quantum system |18j . and therefore highly entangled Floquet 
states are advantageous also in presence of decoherence [29] . 

To quantify the entanglement £j of a Floquet state we consider the time 

average of a predefined entanglement measure |30) £ over one period T: 



(This is reasonable, since is T-periodic quantity.) The particular choice of 

£ will depend on the number of the two-levels under investigation in the subsequent 
sections. Irrespectively of this choice, however, one always has maximal Floquet 
state entanglement £ j = 1 if and only if |$j(t)) is maximally entangled for all times 
t G [0, T), since every (normalized) entanglement measure is a non- negative function 
that vanishes for separable states and takes its maximum value £ = 1 for maximally 
entangled states. Vice versa, Floquet state entanglement vanishes if and only if |$i(f)) 
is separable for all t £ [0, T). 

2.3. Mathematical structure of the Floquet problem 

Before we investigate Floquet state entanglement in detail, we elaborate on the 
mathematical structure of the Floquet problem, which consists in finding the Floquet 
states |$j(t)) and respective quasi-energies £j, for a given time-periodic Hamiltonian 
H(t). This discussion will help us to analyze the phenomena encountered later. 
We start by inserting ^ into Q, what leads to 



This equation is reminiscent of the eigenvalue problem H\&i) = Ei\$i) for 
autonomous quantum systems, and implies that Floquet states and quasi-energies 
are eigenvectors and eigenvalues of the Floquet Hamiltonian Hp = H{t) — idt- 
Differently to the static case, however, Hp acts on the extended Floquet Hilbert space 
Hp = H <S> L 2 ([0, T)) of all T-periodic orbits in the original Hilbert space H. Since 
the space L 2 ([0,T)) of all T-periodic, square-integrable functions is isomorphic to £ 2 
[51] . it is spanned by a discrete basis set, e.g., {e~ lku)t , k € Z}. Hence, every periodic 
orbit !$;(£)) can be expanded in a discrete Fourier series 




(3) 




(4) 




(5) 



with 




(6) 



The periodic orbit is then represented by a "double ket" in Hf, 




(7) 
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where \k) denotes the Fourier basis functions e lkult of L 2 ([0,T)). After expressing 
the Floquet Hamiltonian Hp in the Fourier basis as well, 



[Hi- k + 8 kl kuj}<Ei\k) (l\, 



k.iei 



with 



Eq. Q reads 



Up pi)) = Ei 



(8) 



(9) 



(10) 



Writing out the Fourier indices and Z in matrix form instead, we see that Floquet 
states and eigenvalues of H (t) can be found by diagonalizing the matrix 



(■ 



Ho 



Hi 

m 



Hi 

Yt 



H 2 
Hi 



H H 



V 



(11) 



■J 



This representation of the Floquet problem has a few interesting implications: 
(i) A major advantage of rephrasing the Schrodinger equation as eigenvalue problem 



in Eq. ( 10 ) is that all the concepts developed for autonomous quantum systems 
can be imported. In particular, time-independent perturbation theory can be 



applied to Eq. (10) [IT]. For example, assume that the Floquet states ||$i)) of a 
certain T-periodic H (t) are known. The impact of a small perturbation V(t) (of 
the same periodicity) is then determined, to first order, by matrix elements 



(12) 



k,l 



in Floquet Hilbert space. In the time domain, this reads 



T 



dt(Mt)\V(t)\^(t)). 



(13) 



In particular, if two unperturbed quasi-energies £; and £j of H(t) cross (as some 
system parameter is varied), V(t) can lift this degeneracy and lead to an avoided 
crossing with minimal level separation 2|cjj|. 



(ii) The diagonal of ( 11 ) contains the static part Ho of the system Hamiltonian, while 



the off-diagonal entries Hj reflect driving with frequency jcu. E.g., in case of an 
autonomous Hamiltonian H with eigenstates |$i) and eigenvalues we have 
Ho = H and all other Hj = 0. Then, |$i) €5 |0) is an eigenstate of Hp with 
quasi-energy E^; but all |$i) €5 |fc) are eigenstates as well, with quasi-energies 
Ei + koj. This reflects precisely the previously discussed rest class structure of 
Floquet states. 
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(iii) Despite the fact that the driving field is not quantized in the Floquet picture, but 
is rather represented by a time-dependent classical field in the Hamiltonian, the 
Fourier index k in Eqs. ([5|-([9]) can be interpreted as the number of quanta w in the 
driving field. In fact, an exact correspondence between the Floquet Hamiltonian 
in the Fourier basis and the "dressed" Hamiltonian of the quantized case can be 
established for large occupation numbers of the quantized mode [ITJ [32] . 

(iv) Being interested in N two-level systems in this paper, the extended Floquet 
Hilbert space is U F = [C 2 ]® N <g) L 2 ([0,T)) throughout the sequel of this paper. 
One can formally interpret this as the Hilbert space of an (7V+l)-partite quantum 
system (the (N + l)-th particle having not just two, but countably infinitely many 
levels), with dynamics generated by the autonomous Hamiltonian Hp. 

(v) The number of non- vanishing Fourier components |$j(fc)) is determined by the 



structure of the driving terms Hj in (11) (see, e.g., Ref. [13]). and decays with 



large \k\. Hence, the Fourier index k can effectively be confined to finite intervals 
—M < k < M. Then, H f becomes a finite (2M2 N ) x (2M 2 N ) matrix, which 
can be diagonalized by standard numerical methods. In |Appcndix A[ we derive 
the rule of thumb M ss 2F/uj for a qubit that is driven by a monochromatic field 
of amplitude F and frequency to. 

3. Two weakly coupled qubits under external driving 

As starting point of our investigation, we consider two two-level systems (qubits) that 
are driven by an external field with periodic amplitude f(t) = /(i+27r/w), and coupled 
by a <7 + ct_ "excitation exchange" interaction. The Hamiltonian reads 

2 

H(t) = £ + f(t)4 n A + H m , (14) 



n=l 



ff qq = c(4V 2) +^M 2) J, as) 

with Pauli operators ct*-™- 1 acting on the n-th qubit, u; the single qubit energy splitting, 
and C the qubit-qubit interaction strength. The term f(t)a^ describes a coherent, 
classical electro-magnetic driving field. We assume it to be identical for both qubits, 
so that no individual addressability of the qubits is required here. 



Hamiltonian ( 14 ) is encountered in various physical scenarios, prominent 
examples of which are listed in Table [l] In case of trapped ions 33J and 
superconducting qubits |34j . entanglement is routinely being measured in experiments, 
e.g., by state tomography. The entanglement of Floquet states, as investigated 
presently, can be demonstrated experimentally in the very same fashion. Note that 
in some of the scenarios of Table [TJ the qubit-qubit interaction H qq is not described 
by a er+<7_ excitation exchange mechanism. However, as long as qubits do interact 
pairwise in some form, the details of H qq are not crucial for the phenomena discussed 



in the following, as we will argue in Section 3.4 



Since the Hamiltonian ( 14 1 does not distinguish one or the other qubit, its Floquet 
states can be chosen as either symmetric or antisymmetric under exchange of the 
qubits. The antisymmetric subspace of two qubits, however, consists of the singlet 
(Itl) _ lit)) on ly- (We adopt the spin-1/2 notation |f) and ||) for the qubit basis 
throughout this paper.) Therefore, this maximally entangled state is always one of the 
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driving 


microwave 


laser or 


microwave 


laser 


source 




microwave 






qubit-qubit 


inductively / 


phonon- 


dipole-dipole 


dipole-dipole 


interaction 


capacitively, 
or via cavity 


mediated 






lu /2tt 


5 - 20 GHz 


1 - 10 GHz 


0.5 - 3 GHz 


2.4-10 14 Hz 


w/2tt 


MHz - GHz 


MHz - GHz 


MHz - GHz 


10 14 - 10 15 Hz 


F/2tt 


50 - 500 MHz 


30 MHz 


30 - 400 MHz 


0.5 - 5 MHz 


C/2n 


80 MHz 


30 Hz 


40 kHz 


50 MHz 



Table 1. Energy scales for various physical realizations of the Hamiltonian | |14| l. 
All values are exemplary and meant to indicate orders of magnitude. (In case of 
trapped ions, numbers refer to microwave driving.) Subject of this article is the 
optimal choice of the amplitude F and of the frequency w of the driving field /(t), 
in order to maximize Floquet state entanglement. 



four Floquet states of the Hamiltonian (14 1, independently of the system parameters 



Consequently, we will restrict our discussion to the remaining three Floquet states in 
the symmetric subspace, which is spanned by {|tt) j ^75 (Itl) + lit)) 1 144)}- 

3.1. Monochromatic driving 

An easily implementable driving scheme is monochromatic driving with frequency ui 



and amplitude F, i.e., f(t) = F cos(wt) in Eq. (14 1. 

Figure [l] shows the entanglement Si (indicated by the color code) of the three 
symmetric Floquet states, determined by numerically diagonalizing the Floquet 
problem ( |11[ ), for different values of the system parameters. As entanglement measure 
£ in definition we choose the concurrence [H] 

W)) = |(^> s ®o- s |V)|. (16) 

The axes of each plot span the parameter plane of driving strength F and qubit energy 
ujq; the three rows show different coupling strengths C. All parameters are measured 
in units to the fourth parameter lu. In this way, we scan the entire parameter range 
of Hamiltonian (14) in Figure [l] 



For weak interaction between the qubits, C = 0.02 u (top row), the third Floquet 
state (right column) is almost always maximally entangled. On the contrary, the other 
two states have vanishing entanglement almost everywhere, except for small values of 
Wo and along narrow ridges in the -F-Wo-plane. As C increases (middle and bottom 
row), these entanglement resonances become broader and begin over overlap. 

In order to understand the rich patterns observed on Figure [T] we start with the 
weakly interacting case (C <ti f(t), ujq), in which the qubit-qubit interaction is a small 
perturbation to the non-interacting scenario of C — 0. The unperturbed Floquet 
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Figure 1. Floquet state entanglement Si of two coupled qubits under 
monochromatic driving. Each plot shows the parameter plane spanned by driving 
amplitude F and qubit energy splitting ojoi both measured in units of the driving 
frequency ui. The entanglement is visualized by a color code ranging from white 
(separable) to black (maximally entangled). The three columns correspond to the 
three Floquet states of the symmetric subspace. (The remaining forth Floquet 
state is anti-symmetric and always maximally entangled.) In the upper row, 
where the qubit-qubit coupling C = 0.02 cj is weak, the first two Floquet states 
have close-to-zero entanglement in large parts of the parameter space, and only 
along sharp ridges do entanglement "resonances" emerge. As C increases in the 
middle and bottom row, the resonances become broader and overlap. 



states in the symmetric subspace read 

!$!(*)) = |«Mt)> ® |«Mt)> , (17) 
|* 2 (t)> = |0_(t)>®|0_(t)), and 

l*a(t)> = 4(I<M*)> ® IM*)> + I'M*)) ® l<M*))), 

where \4>±(t)) denote the two Floquet eigenstates of the single qubit Hamiltonian 

h(t) = + Fcos(ujt)a x . (18) 

|$i(f)) and I *&2 (*)} are separable at all times, hence £ \ = £ 2 = 0. |<&3(f)), on the other 
hand, is local unitarily equivalent to -L(|tl) + lit)); and hence maximally entangled 
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at all times, implying £3 = 1. This local unitary equivalence is established by a 
time-periodic transformation U^t) ® U^it), with 

^(*) = It) <0+(*)l + 1+) <<M*)I ; (is) 

the unitarity of U<p(t) is guaranteed by the fact that the two single qubit Floquet states 
\4>±(t)) are orthonormal at all t. To summarize: At C = 0, the left and central plot 
of Figure [T] would appear entirely white, and the right one entirely black. 

As long as the unperturbed Floquet states have non-degenerate quasi-energies, 



perturbation theory (applied to the Floquet Hamiltonian H^?, as discussed in Sec. 2.3| 
guarantees that their character is not drastically altered in the presence of a weak 
qubit-qubit interaction C. I.e., their entanglement should approach smoothly the 
non-interacting values in the limit C — > 0. This reasoning explains why E\ and £2 
vanish in large parts of the parameter plane in the upper row of Figure [l] and why £3 
is predominantly maximal. Only in the vicinity of degeneracies of Floquet eigenvalues 
is a deviation from this picture possible. The resonant behavior of entanglement 
observed in Figure [T] must therefore be anchored in near-degeneracies of some quasi- 
energies. 

Thus, in order to understand the shape and position of the entanglement 
resonances in Figure [l] we have to study the Floquet spectrum of the non- interacting 
system. To this end, it suffices to know the two quasi-energies /i + and /i_ of the single 
qubit Floquet problem 

[/i(t)-i£> t ]|0±(<))=M±l<M*))- (20) 

Since the sum of quasi-energies always equals the trace of the time-averaged 
Hamiltonian [TT] . 

/x+ + At- = — / Tr[h(t)]dt (mod w), (21) 
27r Jo 

and Tr[/i(i)] = in our case, the two single qubit quasi-energies /i± are not 
independent, but fulfill fi + = —fi- (if we work in the central Floquet zone {—%■, 
Therefore, we will speak of the single qubit quasi-energy fi = fi + = —fi- in the 
following . 

With this, the quasi-energies £j of the unperturbed two-qubit Floquet states 



\§i{t)) of Eq. (171 read e\ = 2/i, £ 2 = —2^, and £3 = 0. Thus, the unperturbed 
levels are degenerate at ^ = 0. But this is not the only possible degeneracy between 
unperturbed Floquet states: If, e.g., [i = w/4, we have 6\ =62+^, and thus \&i{t)) is 
degenerate with the shifted Floquet state e lut | )) , since the latter has quasi-energy 
£2 + In general, degeneracy occurs whenever e\ = £2 + nuj, i.e., at 

Ai = n^ (n G No). (22) 

Only in regions of the parameter space where this condition is (approximately) fulfilled 
will the Floquet states significantly deviate from those of the unperturbed, non- 
interacting system at C = 0. 

Figure [2]j a) visualizes condition (22), by quantifying the deviation of \i from the 



nearest n-photon resonance condition, in the F-wp-plane; \x is obtained by numerically 



solving Eq. (20). As expected, the patterns in this plot reproduce the shape of the 
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(a) (b) 




1 2 3 1 2 3 

Driving strength F / eo Driving strength F / co 



Figure 2. Detailed analysis of the "entanglement resonances" observed in 
Figure [l] for weak qubit-qubit coupling C = 0.02 a;. (a) Deviation from the 
degeneracy condition |22| , quantified by min n |4/z/u — n|, as a function of driving 
strength F and of the qubit energy splitting uig ■ The color code ranges from white 
(exact degeneracy) to black (quasi-encrgics as distant as possible), (b) £\, the 
entanglement of the least entangled Floquet state in Figure [l] superimposed by 
contours (dashed blue lines) of (a). Within these contour lines, the degeneracy 
condition is fulfilled up to a finite tolerance in the order of C. This procedure 
predicts the position of the resonances accurately, except for the fact that one out 
of two "corridors" enclosed by the contour lines contains no resonance. 

entanglement resonances. This is verified by Figure [2jb) , which shows the same data 
as in the top left panel of Figure [l] superimposed by contour lines extracted from 
Figure [2ja) . The contour values are chosen such that the lines form corridors in 
which the deviation from degeneracy is small (i.e., less than the interaction strength 
C). Thereby, we obtain an accurate description of the position and shape of the 
resonances. Note, however, that the approximate degeneracy of quasi-energies is only 
a necessary, but not a sufficient criterion for entanglement resonance, since one out 
of two corridors in Figure [2jb) contains no resonance. The absence of resonances at 
these degeneracies will be discussed in Sec. |3.2| 

Figure [3] analyzes the situation in more detail. It shows a section through the 
parameter plane at fixed qubit energy ujq = uj. Both the Floquet spectra of the 
weakly interacting system at C = 0.02 uj (lower panels, solid black lines), and of 
the non-interacting system at C = (dashed red lines) are plotted, along with the 
entanglement £\ (upper panels). As expected, entanglement peaks in the vicinity 
of near-degeneracies of the spectrum. The magnification of a crossing region (right 
panels) reveals that the interacting levels E\ and £2 avoid to cross, with a minimal 
energy distance on the scale of the interaction strength C = 0.02 uj. Such avoided 
crossings of quasi-energies are, in fact, the key to explain the resonant behavior of 
entanglement: Far away from their center, the Floquet states |$i(i)) and |$2(*)) of 
the weakly interacting system are well described by the separable "unperturbed" states 
\<t>±(t))®\<f>±(t)) . At the center of the avoided crossing, however, they become balanced 
superpositions of these product states (in complete analogy to avoided crossings in 
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Figure 3. Section through the parameter space of Figure [2] at fixed qubit level 
splitting ojo = d. Bottom left panel: two Floquet zones of the quasi-energy 
spectrum; both the levels of the weakly interacting system (black solid lines) and 
their non-interacting approximations (red dashed lines) are plotted. The two sets 
are virtually identical on this scale, due to the weak interaction C = 0.02 u. The 
panel above verifies that Floquet state entanglement Si peaks precisely in vicinity 
of (near) degeneracies of the spectrum. Right panels: Magnification of a resonant 
region, revealing an avoided level crossing in the quasi-energy spectrum, on a 
scale in the order of C (bottom, only levels of the interacting system are plotted). 
£l peaks precisely at its center (top). The character of the repelling levels is 
indicated by a label in spin-1/2 notation; e.g., as explained i n th e text, the upper 
level continuously changes its character from 1 3*2 (*))•. c£. Eq. (17) , which is locally 
unitarily equivalent to |44-)i to |<E>i(i)}, which is equivalent to Iff). Note that the 
flat level £3 is not involved in the avoided crossing scenario, since the matrix 
elements | C13 1 and | C23 1 describing its coupling to the other states, cf. Eq. (26) , 
are vanishing here. 



spectra of autonomous quantum systems [45]): 

l*i(*)) - 4(l^+(*)> ® I'M*)) + IM*)) ® |<M*)», 
|* 2 (t)) = -^(\Mt)) ® l<M*)> - \<P-(t)) ® l-M*)))- 

These are maximally entangled states at all times t, as can be seen, once again, by 
application of the local unitary transformation U$(t) ® U$(t), Eq. (19 1. 

Precisely the same mechanism occurs at higher order resonances, i.e., when 
condition ([22]) is fulfilled for n > 0. Then, |$i(t)) is degenerate with e mu3t |0 2 (*)) - 

3etw 

(|0+(t)) ® |0+(t)) ± e iwt \<P-(t)) ® |^_(t))) 



Hence, at the center of an avoided crossing between these levels, we have 
1 

71 



as Floquet states of the interacting system, which again are maximally entangled 
states at all times 

From our perturbative analysis, one can also understand what happens in the 
case of larger qubit-qubit coupling strength: As C grows, the avoided crossings widen 

|| To see this, one has to modify Uj,(t) of Eq. ^ to Cfy(t) = |t> (<£+(*) | + e" 171 "*/ 2 \i) (<j>-(t) |. 
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and migrate from their original position, resulting in broader and slightly shifted 
entanglement resonances that eventually overlap, and form the rich patterns observed 
in the middle and bottom row of Figure [T] 

To summarize our discussion so far, we have identified the central mechanism 
the underlies entanglement resonances: Whenever the weak interaction between the 
qubits lifts a degeneracy of the non-interacting two-qubit Floquet spectrum, it locally 
strongly couples the anti-crossing Floquet states and transforms them from separable 
into maximally entangled states. For clarity, let us emphasize that this mechanism is 
not unique to avoided crossings in Floquet spectra, but similarly occurs for eigenstates 
of weakly interacting, autonomous quantum systems, whenever the corresponding 
energy levels cross under variation of a static control parameter. This phenomenon 
has been discussed, e.g., for spins chains [22l [23] . However, control of entanglement 
(and more generally, of TV-body interaction) by oscillating fields, as suggested here, 
is much more versatile, since the rest class structure of the Floquet spectrum leads 
to multi-photon resonance conditions (alike (22 1), and thereby allows to address the 
quantum many-particle system through a multitude of side-bands, which are absent 
in static control scenarios, and which might experimentally be much easier to access. 



3.2. Detailed analysis of entanglement resonances 

Based on our above, qualitative explanation of the entanglement resonances, we now 
develop a better understanding of their position and line shape in parameter space. 
This involves two steps: 

(i) First, we derive (approximate) expressions for the single qubit quasi-energy \x 
as a function of the driving amplitude F and of the qubit energy splitting ujq. 



This allows to parametrize the degeneracy condition (22), and thus provides an 
analytical description of the contour lines of Figure [2](b) . 

(ii) Next, we study the coupling strength between Floquet levels, i.e., the width of 
avoided crossings which generate the entanglement resonances, and investigate 
why it vanishes within those corridors that contain no resonance. 

The first step can be achieved only within certain approximations, since the single 



driven qubit problem, Eq. ( 20 ) , has no closed solution for the monochromatic driving 
scheme of Eq. (18). In the regime of weak driving, F < tuo, and small detuning, 
ojq ~ lu, the rotating wave approximation (RWA) can be applied [46]. Note that these 
conditions are fulfilled only in a small fraction of the parameter plane shown in Figs. 
[I] and [2j namely in the vicinity of (F = 0,ujq = lo) = Pq. The RWA consists of 
neglecting the "counter-rotating" terms of the driving field, i.e., in replacing 



Fcos(u;t)a x = ~ (e lut a+ + e~ iut a + + e Mt a_ + e~ iut uJ) 



by 



F 
2" 



e <7-l + e <r_ 



This way, the Floquet Hamiltonian ( |TT] ) decomposes into decoupled 2x2 blocks (each 
spread out over two Fourier components) and can be diagonalized by hand, leading to 
quasi-energies 



,,RWA 



±\ (w + V(^-^o) 2 + ^ 2 ) • (23) 
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With this, the degeneracy condition (22) turns into 



yj(u- wo) a + F 2 = nu/2 (neN ). (24) 

This defines a circles of radius nw/2 around Po in the P-wo-plane. Hence, for n = 0, 
condition (24) predicts a point-like resonance at Po- It is hardly visible in Figure [2](b) , 
but appears broader in the middle row of Figure [T] For n = 1, the condition describes 
a circular corridor of radius uj/2 around Po, as visible in Figure [2jb). Beyond this 
corridor, the RWA is no longer justified, and, accordingly, the location of all other 
corridors in Figure [2jb) deviates quite strongly from the circular shape. Only in 
the weak driving limit of higher order resonances, i.e., at P <C u> and ujq ~ 2lu, 
ujq 3w, etc., does RWA provide a reasonably accurate description. In these 
parameter regions, multi-photon transitions are resonant with the unperturbed qubit 
transition frequencies, and the Floquet Hamiltonian again effectively decouples into 
2x2 blocks, with the same degeneracy condition ( p4| [TTj . 

For large driving strength P, the Floquet Hamiltonian strongly couples more than 
two frequency components, and can no longer be approximated by decoupled 2x2 



blocks. This leads to a failure of the RWA resonance condition (24). The Floquet 
states acquire a large number of non-vanishing Fourier components, and this makes 
it difficult to analytically explain the line shapes of the resonances. However, at least 
in the regime of small qubit energies ujq <C w, an approximation can be found that 
is valid for arbitrary driving strengths |llj : It is given in terms of the zeroth order 
Bessel function of the first kind, J : 

^ «. = ^o Jo ( 2 l\ _ (25) 



2 V ui 



This expression vanishes at u>o = 0, consistent with the degeneracy observed in Figure[2] 
in the limit ujq — > . Furthermore, it accurately describes the resonance positions close 
to the P-axis, which coincide with the zeros of Jq. In Figure |4j the exact contour lines 



are plotted along with the approximate expressions (24) and (|25j), to illustrate their 
respective regions of validity. 

The widths of entanglement resonances are given by the widths of the associated 



avoided crossings. According to Eq. (13), this width is determined, at first order 



perturbation theory, by the matrix elements c^- of the participating Floquet states 
|3>i(t)) and \$j(t)) with respect to the qubit-qubit interaction Pf qq : 



2tt 



alt (^(t)\C(a^a^ + a™ a™ . (26) 



Thus, at first order, the width of the resonances scales linearly in the qubit-qubit 
interaction strength C and is bounded by 2|cy| < 2C. This is why the contour lines in 
Figure |2jb) are chosen such that they enclose regions in which the unperturbed levels 
come closer than 2C, thereby defining an upper limit for the width of the resonances. 

The absence of entanglement resonances in some of the such defined corridors, as 
observed in Figure 2|b), can be traced back to a generalized parity of the single qubit 
Floquet problem ( 20j) fTTJ Q21 [5^1 E7] . The operator Sp associated with this parity 
symmetry acts on the Floquet Hilbert space of a single driven qubit and reads 



S P \4>(t)) =<y z \4>{t + n/u)) 



(27) 



Entanglement resonances of driven multi-partite quantum systems 



14 




Figure 4. Illustration of two analytical approximations to the contour lines of 
Figure [2|b). For weak driving amplitude F, a description in terms of the RWA 
holds, explaining the circular shape of contour lines in this area, cf. Eq. \2A\ . For 
small qubit energies u)q <C F, the single qubit quasi-energy fj, can be approximated 
by a Bessel function, cf. Eq. j25| l. The two approximations (bold black lines) 
render correctly the shape of the resonances in the respective regions of parameter 
space. For large qubit energy and strong driving, no analytical description is 
available, and one has to resort to the numerical determination of fi (thin blue 
lines) . 




Figure 5. Three Floquet zones of the Floquet spectrum of two 
monochromatically driven qubits without qubit-qubit interaction (C = 0), at 
ujq = u). Levels are colored red and blue, according to their symmetry with 
respect to the generalized parity iSp2, Eq. ( |28[ |. Since the symmetry is retained 
in the presence of qubit-qubit interaction, level crossings between different classes 
do not turn into avoided crossings in presence of a finite interaction (C > 0). 
This explains the absence of entanglement resonances in one out of two corridors 
in Figure [2jb) . For clarity, we emphasize that the level distortion around 
F = 1.5...2cj is not related to the interaction-induced avoided crossings, which 
lead to entanglement resonances (after all, the plot is obtained for C = 0). This 
distortion is rather a property of the single driven qubit spectrum for strong 
driving, indicating the departure from the RWA regime. 



Since it fulfills (Sp) 2 = 1 and commutes with the Floquet Hamiltonian h(t) — idt, the 
single qubit Floquet states \(j) + (t)) and \(f)-(t)) are eigenstates of Sp, with eigenvalues 
±1, respectively. In the two qubit Floquet Hilbert space, the corresponding parity 
Sp2 reads 

S Pa |$(f))=<72o2|$(i + 7r/w)). (28) 
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The non-interacting Floquet state |$i(i)) = \4>+(t)} , cf. Eq. (17), has parity +1 
with respect to Sp2, because 

s P2 \Mt)f 2 = (s P \Mt))f 2 = (+i) 2 \Mt)f 2 ■ (29) 

Likewise, the second non-interacting Floquet state | *&2 (^)) = 10- (t))® 2 has positive 
parity, whereas \$3(t)) — -^(\<j) + (t)) <£> \<f>_ (t)) + \cf)-(t)) 65 \(j) + (t))) has negative parity. 
Since also the qubit-qubit interaction term that we considered so far (the ct+ct_ 
"excitation exchange interaction") commutes with Sp2, levels of different parity are 
not coupled by this interaction. Consequently, driving-induced degeneracies of two- 
qubit Floquet states of opposite parity are not lifted by a non-vanishing qubit-qubit 
interaction strength C. Only levels of equal parity are coupled by H qq and can thus 
give rise to an entanglement resonance. 

To illustrate the above discussion, three Floquet zones of the non-interacting 
quasi-energy spectrum are shown in Figure [5j with levels colored according to their 
parity. As explained in Sec. |2.1[ the rest class structure of Floquet states results in 
an w-periodicity of the Floquet spectrum. This implies that for each Floquet state 
|<&i(i)) of quasi-energy e i; there is a homologue Floquet state | e wt of quasi- 

energy £j — oj in the neighboring Floquet zone. If |$i(i)) is symmetric with respect 
to Sp2, then |$j(t)) e lLJt is antisymmetric, and vice versa, since of e tu >( t + 7r 1^) = — e 2 "*; 
i.e., the parity of a level switches from Floquet zone to Floquet zone. This implies 
that an avoided crossing of levels is symmetry-forbidden whenever the degeneracy 
condition ( |22[ ) is fulfilled for odd n, and explains why precisely every other corridor in 
Figure J2jb) contains no entanglement resonance. 

3.3. Bi-chromatic, saw-tooth, and 6 -kicked driving 

To underpin our discussion of symmetry-suppressed resonances, we consider a bi- 
chromatic driving scheme, f(t) = F cos(u))t + F' cos(2ujt), which is easily generated 
in laboratories by using the second harmonic field mode [48 . This scheme breaks 
the generalized symmetry Sp. Therefore, all levels couple to each other via the 
qubit-qubit interaction, and we expect an entanglement resonance whenever the 
degeneracy condition (22) is fulfilled. Figure [6] shows the entanglement £\ of the 



first Floquet state, and confirms that now all corridors defined by the degeneracy 
condition enclose an entanglement resonance. (We omit plotting £ 2 and £3 for the bi- 
chromatic case, since, as in the monochromatic case, £2 is very similar to £\, while £3 
is maximal almost everywhere.) Compared to the monochromatic case, the corridors 
are parametrized differently here, since they derive from the single qubit quasi-energy 
ji under bi-chromatic driving, which is different from the monochromatic case. Still, 
our prediction scheme based on the degeneracy condition accurately describes the 
resonances. 

Next, as an example of analytically solvable single qubit dynamics, we consider 
a saw-tooth driving profile f(t) = F ■ mod 1 — |] that periodically ramps up the 
driving amplitude from —F/2 to +F/2, with period T = 2it/lu. Thus, the single qubit 
Hamiltonian is 



h(t) = fa z 



Tf, mod 1 - - 



(30) 



Note that this describes a periodic repetition of a Landau-Zener scenario (in which the 
diabatic states are the eigenstates of a x , and uj plays the role of a coupling strength 
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(a) F' = F/2 (b) F' = F 




Driving strength F / to Driving strength F / co 



Figure 6. Entanglement E\ of the least entangled symmetric Floquet state 
of two weakly coupled qubits (C = 0.02 uj) under bi-chromatic driving, i.e., 
f(t) = Fcos(cj)t + F' cos(2o)t). (a) F' = F/2. (b) F' = F. The predictions 
(dashed blue lines) are contour lines derived from the degeneracy condition 
as explained in the caption of Figure [2jb) . 



(a) Sawtooth driving 




l 2 

Driving strength F / to 



(b) Driving by periodic <5-kicks 

3| 1 • i • 1 1 ■ 



| 



o 



0.5 1 
Driving strength F / io 



1.5 



Figure 7. Entanglement £± of two weakly coupled qubits, driven by (a) a saw- 
tooth profile f(t) = F ■ [± mod 1- \], (b) periodic 5-kicks, f(t) = FT ■ 8(t 
mod T). The predictions (dashed blue lines) are obtained from the exact 
analytical expressions ( |31| and ( |32| , respectively. 



between the diabatic states). Based the analytical solution of a single Landau-Zener 
transition |49j . we derive the explicit expression 



1 

a = — arccos 
p T 



-iuj$T 1 iFT\ 



16F 



) 



(31) 



for the single qubit quasi-energy /j, in |Appendix B[ with iF% denoting Rummer's 
function |50j . Inserting this into the degeneracy condition ( 22 1 , one accurately predicts 
the positions of entanglement resonances, as illustrated in Figure [7]ja). 

In the same Appendix, we also derive the analytical expression for the quasi- 
energies of a qubit exposed to a train of <5-kicks [51] , i.e., for the driving profile 
f(t) = FT-S{t mod T): 



a = — arccos 
' T 



cos (^pj ■ cos(FT) 



(32) 
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Figure 8. Influence of the qubit-qubit interaction H qq on the Floquet state 
entanglement Si of two weakly coupled, monochromatically driven qubits. The 
left panel shows _ff qq = C(p^a 2 _ + a^al), with C = 0.02 cj, as considered in the 
previous figures (i.e., it is identical to Figure[2fb)). In the middle and right panels, 

the interaction term is replaced by C and C^+aFjl^'+^VZ, 

respectively. The latter interaction breaks the generalized symmetry <Sp2, and 
therefore leads to twice as many resonances. 



Also here, entanglement resonances appear in the parameter regions predicted by the 
degeneracy condition, cf. Figure ^o)- 

Generally speaking, one can create entanglement resonances with virtually any 
driving profile f(t), if only the single qubit quasi-energy \i can be tuned to zero by 
variation of the driving parameters. 

3.4- Variation of the interaction term 



So far, we studied Hamiltonian (14) for different driving profiles /(£), but fixed 
the interaction mechanism H qq between the qubits. Based on our understanding of 
entanglement resonances, we expect the details of this interaction not to be decisive for 
the phenomenon: the corridors that define parameter regions of driving-induced near- 
degeneracies of the uncoupled Floquet states - i.e., the regions where entanglement 
resonances may arise - are independent of H qq , since they derive from the degeneracy 
condition (22), which relies on the single qubit quasi-energy fi alone. The role of H qq 



is to lift the degeneracies of the non-interacting spectrum, and any generic two-qubit 
operator will achieve this. Of course, the details of the interaction determine the exact 



value of the coupling matrix elements in Eq. (26), and, hence, define the width 
of the resonances; furthermore, we have seen above that, if the interaction preserves 
a certain symmetry of the single qubit Floquet Hamiltonian, some degeneracies may 
not be lifted, i.e., some entanglement resonances may be symmetry- forbidden. 

This reasoning is illustrated in Figure [8] It shows Floquet state entanglement 
under monochromatic driving for different interaction mechanisms. The left plot 



corresponds to the excitation exchange interaction of Eq. (15). It is replaced by 
ai 1} ai 2) in the central, and by (aP + ai 1} ) (c4 2) + cri 2) )/2 in the right plot [fj The 

1 If the Pauli matrices refer to a physical spin-1/2 particle, and not to an abstract two-level system, 

the different interactions have the following physical meaning: cr^a^ and (cr^ +ct# ')(&^ +<?^' ') /2 
both describe dipole-dipole interactions between the spins. In the former case, the axis connecting 
the two spins is perpendicular to the magnetic field leading to the Zeeman term ^po" z in \14\ \ 
in the latter case, the axis connecting the spins and the magnetic field are aligned under an 
angle of 45°. The excitation exchange interaction cr^'cr^ 2 ' + a^a^\ on the other hand, is an 

approximation to <r^a^ that neglects the magnetization-changing terms ax', since <r^a^ = 
(1) (2) , (1) (2) , (1) (2) . (1) (2) 
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left and the central plot of Figure [8] are virtually identical. (Only the resonance 
widths are slightly broader in the central plot.) In particular, the symmetry-forbidden 

resonances remain suppressed, as ci^ci 2 ' commutes, just like the excitation exchange 
interaction, with the generalized parity Sp2- This is no longer the case in the right 
figure, where the interaction breaks the symmetry, thereby coupling all levels to each 
other, and triggering entanglement resonances at all driving-induced degeneracies of 
the uncoupled Floquet spectrum. 

To summarize, we have seen that the observed entanglement resonances are largely 
independent of the interaction mechanism. They can accurately be predicted by 



solving the single qubit Floquet problem, Eq. ( 20 ) , alone. This is by far less demanding 
than solving the full Floquet problem for two qubits - a fact that becomes even more 
advantageous in the case of more than two qubits, that we consider in the following. 

4. Three Qubits and GHZ entanglement 

Based on the understanding of Floquet state entanglement of two periodically driven 
qubits, we extend our investigation to three qubits in the following, to see whether also 
in this case entanglement behaves strongly resonantly in certain parameter regions, 
and to pave the way for studying larger numbers of qubits. 

While for two qubits, the definition of a maximally entangled state is unique, 
this is no longer the case for N > 2 qubits [52] • For three qubits, there are two 
inequivalcnt classes of maximally entangled states [S3]: GHZ-entangled states, which 
can be transformed into 

|GHZ) = i (|ttt) + \U\)) (33) 

by local unitary operations, and W-entangled states, that are locally unitarily 
equivalent to 

iw) = ^(iut> + iw> + im))- (34) 

Vo 

By fixing the entanglement measure, one specifies the type of entanglement that is 
quantified. Here, we focus on GHZ entanglement, and use the three-tangle [M] as 
entanglement measure in Eq. |3]), which is maximal for GHZ states and vanishes not 
only for separable and bi-separable state^J but also for W states. 

4-1. Identical system parameters 



A natural extension of the two qubit Hamiltonian ( 14 1 to three qubits is 

3 



H(t) = ]T (^crW + f{t)*M) + H qq , (35) 



71 = 1 

3 



#qq = C (4" )ctM + 4" )(jM ) • (36) 

n<m 

~*~ Bi-separable states are not fully separable, but only separable with respect to a certain bipartition. 
An example is ^=(|ttt> + |Ut» = ^(Itt) + |W» ® It}- 
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Here, the interaction _ff qq couples only two qubits at a time. This is a reasonable model, 
since three qubit interaction is typically much weaker than pairwise interaction (e.g., 
in resonant excitation exchange between Rydberg atoms [55] or chromophores [55]). 
or it is the only relevant interaction by construction (think of engineered systems, 
such as inductively coupled superconducting qubits [57]). On the other hand, the 
assumption of all qubits to have identical coupling strength C and energy splitting ujq 
is somewhat artificial and often violated in reality; we discuss perturbations of this 
idealized situation in the next section. 

For the time being, however, H(t) is invariant under cyclic permutation of the 
qubits. We exploit this three-fold symmetry to divide the Floquet states into three 
decoupled symmetry classes, similar in spirit to the two qubit case, where we split off 
the one-dimensional antisymmetric subspace. Again, we take a perturbative approach 
in the qubit-qubit interaction, and find the symmetric Floquet states of three non- 
interacting, driven qubits, analogously to Eq. ( |17[ ), to be 

|<ih(t)) = |<Mi))|<M*))l<M*)), (37) 
|<h(t)> = |0_(t)>|0_(t)>|<Mt)>, 

l*a(t)> = ^(|<M*))|4> + (*))I<M*)) 

+ \<t> + (t))\<t>-(t))\Mt)) 

+ |0_(*)>IM*)>IM*)»> and 

l*4(*)> = i(|<M*))l<M*))l<M*)) 

+ IMW+(*)> IM*)) 
+ |<M*)> l<M*)> l<M*)»- 



Applying the local unitary transformation U^t) defined in Eq. (19), we can read off 
the entanglement properties of these states: 

[C^(*)]® 3 |$i(*))= |ttt) 

[^(*)]® 3 i<h(t)> = \m) 

pm 93 l*a(*)) = ^(itn) + mt) + urt)) = iw> 
[u^(t)f 3 i$ 4 (t)) = j=mt) + \m) + im» = iw) . 

Both |W) and |W') are maximally W-entangled (since the latter is transformed into 
the former by interchanging [f) and |4<) for each qubit), and altogether we have 
two separable and two W-entangled Floquet states in the permutation-symmetric 
subspace. Hence, all Floquet states have zero GHZ entanglement in the non- 
interacting case. However, if two of these states are tuned into resonance by a 
suitable driving field, a non- vanishing particle-particle interaction will induce resonant 
coupling, and the corresponding Floquet states will turn into superpositions of the 
non-interacting states (37), what offers the possibility for GHZ entanglement. 

Besides the permutation-symmetric subspace, there is the subspace associated 
with the permutation eigenvalue a — e+ 27 ™/ 3 , which contains, in the non-interacting 
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Figure 9. Floquet state entanglement of three monochromatically driven qubits, 
as a function of driving amplitude F and qubit energy splitting u>o, for weak 
interaction between the qubits (C = 0.02 u). The Floquet state with the highest 
amount of GHZ entanglement is shown, in contrast to the previous figures for two 
qubits, in which the entanglement £i of the least entangled state is of interest. 
Entanglement resonances appear exactly in the same parameter regions as in the 
two qubit case. 



case, the Floquet states 

[U<p(t)}® 3 l*6(*)> - ^(Ittl) + a* \Ut) + a \m)) and 

lu^t)}® 3 i$ 6 (t)) = i(iwt) + a* \m) + a im», 

and the subspace associated with the eigenvalue a* = e - 27 ™/3 5 with Floquet states 
obtained from those above by exchange of a and a* . Hence, all unperturbed Floquet 
states in the non-symmetric subspaces are W states. Their superposition yields at most 
very poor GHZ entanglement (the three-tangle is bounded by 1/9 in these subspaces^}, 
and therefore we focus our discussion on the symmetric subspaces in the following, 
which exhibits richer phenomena. 

The results for weak interaction strength C and monochromatic driving, f(t) = 
F cos(wi), are shown in Figure [9] Only the permutation-symmetric Floquet state with 
the highest amount of entanglement is shown; i.e., at a given position in the parameter 
plane, we maximize over all Ei. We observe vanishing entanglement in most parts of 
the parameter plane, but resonant behavior in the same areas as in the two qubit 
case. Accordingly, the phenomenology can be explained in the framework that was 
established in Sec. [3] In the absence of interaction, none of the Floquet states exhibits 
any GHZ entanglement, as reasoned above. In the presence of weak interaction, the 
character of Floquet states is not severely altered, unless the quasi-energies of the 
non-interacting system are near-degenerate. For three qubits, the unperturbed quasi- 
energies are £i = 3/z, £2 = — 3/i, £3 = /i, and £4 = Hence, whenever the degeneracy 



condition /i = nw/4 of Eq. (22 1, which was derived for two qubits, is fulfilled, e\ and 
£4 (and likewise £2 and £3) are degenerate. As in the two qubit case, the generalized 
parity suppresses the coupling for odd n; but for even n, the interaction matrix element 

* This statement is derived as follows: Any superposition a\^>^(t)) + /3|<I > e(t)} is locally unitarily 
equivalent to a |W)+/3 |W). (The corresponding transformation, after application of [[/^(t)]® 3 , reads 
e !7r(o- z — ct z )/3 j rpj^ three-tangle can be explicitly evaluated for this superposition, and is ^|a/3| 4 . 
Hence, it takes its maximum 1/9 for a = f3 = ^g. 
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Figure 10. Level structure of three non-interacting, driven qubits at uiq 
For more details, see caption of the respective two qubit plot in Figure [5] 



C14 is finite, and an avoided crossing opens up under finite qubit-qubit interaction. In 
the vicinity of the avoided crossing, the Floquet states turn into superpositions of the 
unperturbed states, a |$i(i)) + (3 |$4(i)), with \a\ 2 + \(3\ 2 = 1. As one sweeps through 
the anti-crossing, |a| and |/3| continuously vary between and 1, and at a certain point 
one always has the particular superposition 

1 |$ 1 (t))±,/f |* 4 (t)) (38) 




Ittt) ± 




= \ul(tT 3 (ittt) ± im) ± im> ± nit)) /2. 

This is an exact GHZ state, what can be explicitly checked by applying the Hadamard 
transformation R = e~ 17Tay / 4 to every qubit. Therefore, Floquet state entanglement 
becomes maximal at the avoided crossing of e\ and £4. Hence, the degeneracy 
condition fi = ncj/4, that describes the position of resonances in the two qubit case, 
also gives the position of GHZ entanglement resonances of three qubits. 

As depicted in Figure |10| there are more level crossings in the unperturbed 
spectrum of three qubits. The reason why they do not cause entanglement resonances 
in Figure [9] is not obvious. E.g., at the center of an avoided crossing of e± and £2, we 
expect the Floquet states to turn into 

1 (!$!(*)> ± 1*2(0)) = Pl(t)f 3 ^= (Ittt) ± HID) , 



which are maximally GHZ entangled. The degeneracy condition for e\ and £2 is 

(n e N ). 



G 



(39) 



(40) 



As shown in Figure 1 1 entanglement resonances do indeed occur in parameter regions 
where condition ( 40 ) is fulfilled for even n (the odd resonances are again symmetry- 
suppressed). But, given the scale of Figure 11 the width of these resonance is 
extraordinarily small. This is why they are not visible on the scale of Figure [9] The 
reason for this is the relevant coupling matrix element, 



C12 = 



2^ 



2tt/u> 



dt (MtT 3 #qqlM*)> 



03 



(41) 
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Figure 11. Magnification of a small area of Figure [9] A narrow entanglement 
resonance appears, triggered by an avoided crossing of ei and 62- The degeneracy 
condition for these levels is given by Eq. i |4()[ ] . Inside the corridor enclosed 
by the dashed blue lines, this condition is fulfilled for n = 2, up to a finite 
tolerance C 2 , reflecting the second order character of the coupling between the 
participating levels. The slight mismatch between predicted and actual position 
of the resonance is due to the fact that the interaction does not only open up 
avoided crossings, but also slightly shifts their position from the crossing of the 
unperturbed levels 



which vanishes exactly for any H qq that mediates at most two-particle interaction, like 



the qubit-qubit interaction considered here, cf. Eq.(36|. (This follows from the Floquet 
state orthogonality (<j> + (t)\(f>-{t)) = 0, which implies ((f) + (t)\® 3 |<M*))® 3 = 

for arbitrary single particle operators X and Y .) Hence, |<!>i(£)) and |<I>2(i)) couple 
only to second order, via the intermediate states |<E>3(i)) and |$4(t)). Therefore, the 
width of their avoided crossing scales as C 2 , which explains the sharpness of the 
resonance in Figure [TT] 

Finally, avoided crossings between £3 and £4 result in Floquet states that are 
locally unitarily equivalent to a superposition of |W) and |W'). As discussed for 



the non-symmetric permutation subspaces above (see footnote on p. 20), such a 
superposition bears only little GHZ entanglement. In addition, the transition matrix 
element C34 vanishes exactly, so that the avoided crossing between these two levels 
opens only in second order. For both these reasons, entanglement resonances between 
£3 and £4 are not detected on the scale of Figure [9] 

4-2. Non-identical system parameters 



The Hamiltonian considered so far, Eq. (35), is invariant under permutation of qubits 



This assumption allowed us to reduce the number of relevant states, by focussing on 
the symmetric subspace. In particular, we assumed that all qubits have equal energy 
splitting cjo, and that the strength C of the qubit-qubit interaction is identical for all 
pairs of qubits. In certain cases (e.g., in engineered quantum systems), this can be 
a valid assumption; but in most naturally occurring situations, one often rather has 
slightly different parameters for each qubit. 

Therefore, we study a slight deviation from identical system parameters in 
Figure [l2j The left panel shows the influence of individual coupling strengths C for 
each pair of qubits. Each term in the qubit-qubit interaction is of different strength 
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1 2 3 1 2 



Driving strength F / CO Driving strength F / en 

Figure 12. Impact of non-identical qubit parameters on the Floquet state 
entanglement of three monochromatically driven, weakly coupled qubits. As in 
Figure [9] only the Floquet state with the highest amount of GHZ entanglement 
is shown. In the left panel, the interaction strength C between different pairs of 
qubits varies in a range of ±10%. This modification does not significantly affect 
entanglement resonances. In the right panel, the energy splitting ojo is different 
for each qubit, again with a spread of ±10%. This affects the resonances in a 
modest, but discernible way. 



now, 

3 

zj r~i \ " / (n) (m) (n) (m)\ 

H m = C anm l°+ - + a + - ) ' 

reflected by the weighting factors a nmi which are chosen as {ot\2 = 0.9, «23 = 
1.0, «3i = 1.1}. The impact of such modification is not apparent, and the 
entanglement resonances are practically not affected (compare Figure [9]). This is 
consistent with our discussion in Sec. |3.4| where we reasoned why a change of the 
interaction term only marginally influences the resonances. The same argument holds 
here: Since the role of H qq is only to open the avoided crossings, a variation of its 
strength affects only the width of resonances, not their shape. 

The situation is different in right panel of Figure [12] Here, the qubit energy 
splittings ujq are individualized by putting 

H(t) = ]T (A. + f(t)4 n) ) + H qq , (42) 

71=1 

with weighting factors {/3i = 0.9, f5% = 1.0, /?3 = 1.1}. Since this modification 
influences the single qubit quasi-energy fi, the position of level crossings in the 
unperturbed spectrum is affected, and, accordingly the shape of the resonances may 
change. Yet, the two dominant resonances Figure [12] retain their overall shape and 
width, and only the smaller resonances split into "twin arches" . 

In conclusion, we find that a modest deviation from the assumption of identical 
system parameters for each qubit does not significantly alter the phenomenology of 
entanglement resonances. 

5. N > 3 qubits and general discussion of entanglement resonances 

We have seen that the phenomenon of entanglement resonance is neither restricted to 
two qubits, nor to a particular driving profile f(t) or qubit-qubit interaction H qq . In 
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the following, we extend our analysis to the general case of N > 3 weakly interacting, 
driven qubits. 

5.1. Identical system parameters 

In the symmetric case of identical qubit parameters, the ./V-qubit Hamiltonian we want 
to study reads 

N 

H(t) = (iM n) + f®°* n) ) + H ^ ( 43 ) 



As discussed for N = 3 qubits in Sec. 4.1 H qq is typically restricted to two-body 



interactions. As an example, we study again the excitation exchange interaction 

N 



and introduce the collective spin operator 

N 

2 



H qq = cY, (4 n) ^- m) + 4 n) ^- m) ) . ( 44 ) 



1 N 

J = 2 E ° in) ^ (n) = {4 n) ,4 n) '^ n) })- (45) 

n=l 

With J± = (J K ± iJ y )/2, the Hamiltonian rewrites 

H(t) = u Q J z + 2/(i) J x + ^(J+J- + J-J+ - N). (46) 

Since H(t) commutes with J 2 , the dynamics conserves the total spin, with values 
N/2, N/2-l, . . . , 0, for even AT, or N/2, N/2-1, ...,1/2, for odd N. The subspace with 
maximal total spin has dimension N + 1, and contains all states that are symmetric 
under cyclic permutation of the qubits. As in the previous chapters for two and 
three qubits, we focus on this symmetric subspace in our discussion of Floquet state 
entanglement. This is not a severe limitation, since it is reasonable to assume that 
an interacting system of identical qubits is initially in a symmetric state, e.g., the 
de-excited state 14.)®^, and can hence only explore the symmetric subspace under the 



dynamics induced by Hamiltonian ( 45 1 . 
Introducing the Dicke states [581 159] 



\N, m) 



N 



-1/2 

JT\if N , me{0,...,N}, (47) 



as a basis of the symmetric subspace (m corresponds to the number of "excitations" , 
i.e., to the number of qubits in the spin-up state), the symmetric Floquet states of the 
non-interacting system, C — 0, read 

\Mt)) = lUl(t)] 93 \N,i). (48) 



U^it) is defined via the single qubit Floquet states, see Eq. (19 1. Consequently, the 
entanglement properties of |$i(t)) are equivalent to those of \N, i). The latter are 
studied in detail in [501 GEO E3- One finds that |7V, 0) = \\)® N and \N,N) = \\)® N 
are the only separable states, and that the remaining Dicke states belong to distinct 
SLOCC classes (at least for i < N/2 or i < (N + l)/2, respectively, since \N,i) and 
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\N,N — i) are connected by the local operation that exchanges f and 4- labels, and 
therefore have equivalent entanglement properties). All these distinct classes have W 
character [BTJj . 

In the presence of a weak qubit-qubit interaction, the Floquet states |4j(t)) are 
not altered significantly, as long the quasi-energies £j of the non-interacting case are 
far from degeneracy. On the other hand, when two levels Si and Sj cross under 
variation of some parameter (e.g., the driving amplitude F or the driving frequency 
oj, which are usually easily tunable), the interaction can lift this degeneracy and induce 
an avoided crossing, under conditions discussed below. By sweeping the parameters 
through such an avoided crossing, any superposition a|$j(t)) + (3 \ $>j(t)) will become 
a Floquet state of H(t) at some point. This way, Floquet states can be "designed" 
to be locally unitarily equivalent to any superposition of two Dicke states \N,i) and 
\N,j), simply by tuning the respective quasi-energies £i and £j into resonance. This 
opens up the possibility of achieving Floquet states entanglement of many more 
classes, compared to the W classes the bare Dicke state belong to. E.g., the balanced 
superposition of \N, 0) and \N, N) results in the ./V-qubit GHZ state. For four qubits, 

| |4, 0) + y§ 1 4, 3) bears "T-entanglement" , which is SLOCC inequivalent to "single- 
excitation W entanglement" (found in |4, 1) and |4,3)), "two-excitation W" (found in 
|4,2)), and GHZ entanglement [52] . 

In order to design Floquet states in this manner, one needs to determine the 
driving parameters that tune Ei and Sj into resonance. Since the non-interacting quasi- 
energies are = [(N — 2i)jJ\ mod lj, with the single qubit quasi-energy /i defined in 
Eq. (20), the resonance condition reads 

t* = n nr" v neN . (49) 
2 0- 3) 

In order to find the right parameters for the desired level crossing, if therefore suffices 
to solve the Floquet problem of the single qubit Hamiltonian h(t), and to identify 
parameters which fulfill the resonance condition. This is by far simpler than solving 
the Floquet problem of the full, interacting ./V-qubit Hamiltonian H(t). 

It remains to determine how strongly the non-interacting Floquet levels are 
coupled by the qubit-qubit interaction: If two levels are not coupled at all, no avoided 
crossing occurs between them, and no superposition of Dicke states can be created in 
the fashion sketched above. Furthermore, if levels couple very weakly, their avoided 
crossing is small, and the driving parameters have to be tuned very precisely to 
establish the desired superposition, a task that becomes experimentally unfeasible 
below a certain threshold. Analogously to the case of two and three qubits in Eqs. ( 26 1 
and (41 1, the coupling strength is determined in first order by the interaction operator 
Floquet matrix element 

alt (^(t)! JT qq (50) 



2k 



If -f/qq mediates only two-qubit interactions, vanishes for \i — j\ > 2. Therefore, 
|4j(4)) and \$>i+2m(t)) are only coupled at m-th order, and the width of their avoided 
crossing scales like C m . Hence, the ability to superimpose these states by tuning 
into the avoided crossing rapidly decreases with increasing m. The creation of GHZ 
entanglement by superimposing |4 (4)) and |4jv(4)) is particularly affected by this 
issue, since there m = (N + l)/2. For three qubits, this leads to the second order 



character of the entanglement resonance shown in Figure 11 
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5.2. Non-identical system parameters 

The analysis in terms of the collective spin operator J and the symmetric Dicke states 
\N,i) above was possible because all N qubits were assumed to be identical in the 



Hamiltonian (43 1. Now, we drop this assumption and introduce individual parameters 
for each qubit: 



N / In) 
' W 



N 

H qq = E C(n ' m) (°+ + v { + ] °- {m) ) (52) 

n<m 

For three qubits, we have phenomenologically studied the influence of such 



modifications on Floquet state entanglement in Sec. 4.2 

In general, individual parameters for each qubit break the permutation invariance 
of H(t). Generically, this leads to completely different entanglement properties of 
Floquet states than those discussed above for the symmetric case. This is due to 



the single qubit problem, Eq. ( 20 ) , being different for each qubit now, and leading to 
individual single qubit Floquet states \4>±\t)) and quasi-energies fi^ n \ Then, Floquet 
states of the non-interacting system are product states 

N 

\^s(t)) = (g)\^(t)), (53) 

n=l 

labeled by a string s of N plus or minus signs, and have quasi-energies 

N 



*=$>„m (b) - ( 54 ) 



n=l 



Let us first explain how the case of identical parameters fits into this picture: There, 
due to /i^™) = n, quasi-energies eg and ep are degenerate, if s and p contain an 



equal number of plus signs. The product states (531 are not a clever choice of 
basis in these degenerate subspaces: Since we later want to include a permutation- 
invariant perturbation H qq , the appropriate basis is rather the one that diagonalizes 
the permutation operator. This way, focussing on one particular permutation class 
(e.g., the symmetric one), all ambiguities are eliminated: We have quasi-energies 
eg = (N — i)fj, (where i is the number of plus signs in s), and the symmetrized Floquet 
states of Eq. p8| . 

Such systematic degeneracies are absent in the case of non-identical parameters. 
Generically, the only basis of the unperturbed system are the product states \$>g(t)} 



of Eq. (53). Restricting the qubit-qubit coupling i? qq to two-qubit interaction again, 



s 



the product states |$?(t)) and \$>p(t)} are coupled by H qq (in first order) only if 
and p differ in at most two entries. Then, in the vicinity of the avoided crossing of 
and ep, the two corresponding Floquet states turn into the superposition of |$j(i)) 
and |$p-(£)). Since all N — 2 subsystems corresponding to identical entries of s and 
p can be factored out in this superposition, only bipartite entanglement between the 
remaining two qubits is generated at this entanglement resonance. 

The limitation to bipartite entanglement for the first order entanglement 
resonances can only be overcome if more than two levels meet in an avoided crossing. 
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Apparently, this is the case in Figure[l2] where plenty of GHZ entanglement is present. 
This can be ascribed to the fact that only a slight deviation from identical parameters 
was studied there, which retain some systematic degeneracies of the symmetric case, 
in which multipartite entanglement is the norm rather than an exception. 



6. Summary and conclusion 

In this paper, we considered weakly interacting qubits that are coherently driven by 
an external, time-periodic field. This situation is realized in many experiments with 
quantum-mechanical two-level systems, such as trapped ions [53], superconducting 
qubits [36] , Nitrogen- vacancy centers in diamond [40] , or cold Rydberg atoms [41] . By 
employing the Floquet picture, we are able to analyze parameter ranges beyond the 
rotating wave approximation (RWA), i.e., strong and off- resonant driving. The main 
goal was to study the entanglement of the Floquet states of such systems. 

For two qubits, we first analyzed the case of identical parameters for both 
qubits, and found that two of the three Floquet states of the symmetric subspace 
are only entangled in certain regions of the parameter space. The occurrence of 
these entanglement resonances was studied in more detail in the following: An 
explanation was given for the emergence of the resonances, based on a perturbative 
treatment of the qubit-qubit interaction, and the avoided crossings induced by 



this interaction. This led to the necessary, but not sufficient condition (22) for 
entanglement resonance that nicely describes the shape of resonances in parameter 
space. Based on this understanding, we studied different time-dependencies of 
the driving field (monochromatic, bi-chromatic, saw-tooth), different interaction 
mechanisms, and the influence of individual parameters for both qubits, and found 
that all these variations can be well accounted-for by our theoretical framework. 

For three qubits, we focussed on the GHZ entanglement class and again found 
entanglement resonances in parameter space. An explanation of this observation was 
given along the same lines as in the two qubit case. Finally, we generalized our 
reasoning to N > 3 qubits, for both identical and non-identical system parameters. 

Altogether, the main results of this work are 

(i) the finding that maximally entangled Floquet states exist at avoided crossings in 
the Floquet spectrum of weakly interacting, periodically driven quantum systems. 
While a similar phenomenon has been reported in weakly interacting systems 
under static forcing [53] [22] , the advantage of the here proposed control scheme is 
that alternating control fields often provide the simplest way to control a quantum 
system (think of, e.g., trapped ions [8], cold atoms [9], or color centers in diamonds 
[64]). In addition, control via time-periodic fields is more versatile, since the 
quantum system can be addressed through a multitude of side-bands, and it 
can lead to increased coherence times [55], suggesting that the entanglement of 
Floquet states might be more robust against decoherence than the entanglement 
of eigenstates of an undriven system. Verification of this conjecture is a promising 
perspective for future studies. 

(ii) a simple prediction scheme for the positions of the avoided crossings (and therefore 
of entanglement resonances) in parameter space, based on the degeneracy 



condition ( 22 ) . Since this condition involves only the quasi-energy /i of a single 
driven qubit, it does not increase in complexity as the number of qubits increases, 
and can even be evaluated analytically for certain driving profiles. 
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Notably, the prediction scheme is not limited to the analysis of Floquet state 
entanglement, since entanglement is not the only property that behaves critically 
in the vicinity of avoided level crossings: In fact, as discussed in Sec. [3j 
the qubit-qubit interaction (when regarded as a perturbation to the non- 
interacting, periodically driven system) can change the character of Floquet states 
significantly only if two quasi-energies are in resonance, i.e., close to degeneracy. 
Hence, by tuning the driving parameters in or out of such a resonance, one 
effectively enhances or suppresses the interaction in a controlled way. Identifying 
the resonant driving parameters is therefore essential for the analysis of any 
interaction-related property of a periodically driven quantum system. Moreover, 
by bringing more than two levels into resonance, one can dedicatedly study the 
effect of three-body, four-body, etc., interaction. Therefore, the scheme presented 
here does not only explain the behavior of Floquet state entanglement, but more 
generally provides a recipe for controlling many-body interactions by periodic 
driving fields, e.g., in nuclear spin systems [65] or Rydberg gases |21j . 
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Appendix A. Rule of thumb for the number of Fourier components 



In this Appendix, we derive a rule of thumb for the number M of Fourier components 
that have to be taken into account in order to appropriately describe the Floquet 
states of a single, monochromatically driven qubit. 

In the language of Sec. |2.3[ we seek M G N such that 

(A.l) 



Vfc e Z with > M : || \(f>i(k)) || < 1 

Here, \(f>i(k)) is the fc-th Fourier component of the Floquet states \4>i(t)} of 



h(t) 

The respective Floquet Hamiltonian hp is 



F cos(uit)a x . 



I 



Uqo~ z + 2uj Fa x 
Fa x uj a a z Fa x 
Fa x uj a z - 



2uj 



V 



7 



(A.2) 



(A.3) 



in Fourier representation. Due of the generalized parity symmetry Sp, cf. Eq. ( |27[ h 
hp consists of two uncoupled blocks h^, according to the two parity classes [|] 



In the notation of Sec. 



2.3 



|t) ® \k) and \i) ® \k + 1} with even (odd) k 



the block of positive (negative) parity comprises the basis states 
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Driving strength F / CD 



Figure Al. Number of frequency components M of the Floquet states of a 
monochromatically driven qubit, as defined in Eq. ( | A. 1 [ t . The value of M is 
indicated by a color code, and plotted as a function of driving strength F and 
qubit level splitting wq. At cjq = 0, one has M pb 2F/u>, as discussed in the text. 
For loq ^ 0, M is always smaller than this estimate. 



A. 



=F^o + 2u F 
F ±luq F 
F Two - 2w 



V 



(A.4) 



7 



For wo = 0, the eigenvalues of this matrix are ju> (with j £ Z). The fc-th frequency 
component of the corresponding eigenstates \\4>j}} is given by \<f>j(k)) = J\j_u{F/ui) 

with J m denoting the ra-th order Bessel function of the first kind. Hence, \\<f>j)) 
is centered around the j-th frequency component, and the different eigenstates are 
frequency shifted versions of each other. It therefore suffices to study the number of 
frequency components M of ||^o)) : From numerics, we find that \ J\k\(F/uj)\ < .05 
holds for k <; 2F/uj, and therefore we have M = 2F/u>. 

Figure \KT\ depicts the behavior of M for wo ^ 0: At any driving strength F, the 
value of M for uiq ^ is always smaller than for ojq = 0. Therefore, the expression 
M = 2F/ui 1 which is strictly valid for luq — 0, always defines an upper bound of M, 
and can hence be used as a general rule of thumb. 



Appendix B. Quasi-energies for saw-tooth and J-kicked driving 

In this Appendix, we derive the exact quasi-energies of a single qubit, driven with 
a saw-tooth profile, or by periodic 6-kicks. To the best of our knowledge, both 
expressions have not been derived in the literature so far. 
In Sec. 



3.3 



we consider the saw-tooth driving profile f(t) = F ■ [A mod 1 — ¥] . 
It describes an external field that is periodically ramped up from —F/2 to F/2, with 
period T = 2tt/uj. Let us consider, for the moment, an additional constant "offset" 
F/2, such that f(t) = F ■ [A mod l] . (I.e., 
to F now.) Within the time interval [0, T), 
Hamiltonian 

CJ 



the driving field is ramped up from 
the evolution is thus governed by the 



h(t) 



Ft 



T 



ex. 



(B.l) 
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Comparing this to the Hamiltonian 



at 



cr z + Va x 



(B.2) 



of a Landau-Zener scenario [59], we see that h(t) and /ilz(^) are basically identical; 
only the roles of a z and (T^ are interchanged. Indeed, by identifying V = uj/2 and 
a = — 2T/T, and defining the Hadamard transformation R — e~ I7rcr «/ 4 , we have 



h(t) = Rh hz (t)R~ 



(B.3) 



tf(0,t) 

a(t) 

6(i) 



(B.4) 



Since R is static and unitary, quasi-energies are invariant under this transformation. 
To derive the quasi-energies of /ilz (t) , we make use of the fact that its time evolution 
operator U(0,t) is explicitly known [15] : 

'a(t) -&(*)' 
v 6(t)* a(t)* 

'zV^ 1 

2aT' 2' ~~2 

1 iV 2 3 -ta 

2 + 2aT' 2' 

with iFi denoting Kummer's function [50 . The eigenvalues of this operator are 



= i-Fi 



-iat 



-iat 2 /4 



, !ft denoting the real part; hence, its eigenphases are ±arccos(9?[a(i)]). 



D ±iarccos(3?[a(i)]) 

The quasi-energies /x± of hhz(t) - and therefore also of h(t) - are the eigenphases of 
U(0,T), divided by T 67}. Inserting the definition of a(t) and the relations V = uj/2 
and a = —2F/T, we finally have 



, 1 

(i± = ±— arccos 



K[i-Fi 



16F 



iFT) e FT ' 2 ) 



(B.5) 



To derive the quasi-energies in the absence of the additional offset introduced 
above, i.e., for f(t) = F ■ [4 mod 1 — |] , we simply shift the time interval [0, T) by 
half a period, since the time evolution generated by the Landau-Zener Hamiltonian 
hjjiit) during the time interval [— T/2, T/2) precisely described the dynamics of saw- 
tooth driving without the offset. The time evolution operator now reads 

[/(-T/2, T/2) = U\0, -T/2) {7(0, T/2) 

\b\ 2 2a*b 

H 2 -\b\\ 

where we used the short-hand notation a and b instead of a(T/2) and b{T/2). From 
the eigenphases ±arccos(|a| 2 — |6| 2 ) of this operator, we find, together with the 
normalization condition \a\ 2 + \b\ 2 = 1: 



-2ab* 



(B.| 



fi± = ± ^arccos (2|a| 



1 



(B.7) 



± ^ arccos 2 iFi 



-iuilT 
16T ' 



iFT\ 



- 1 



Likewise, the quasi-energies of a periodically (5-kicked qubit can be obtained 
analytically. The driving profile is f(t) = FT ■ 5(t mod T), and the time evolution 
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operator reads 

a(t) = e - *^* cos(FT) 
b(t) =ie^*sin(FT), 
leading to quasi-energies 

/i± = ± ^arccos(5i[a(T)]) 

= ± — arccos 
T 
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